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V  This  report  deals  with  two  problems  in  the  area  of  stochastic 
systems. 1  First,  the  eigenvalue  problem  which  is  solved  using  some  tech¬ 
niques  of  functional  analysis.  Secondly,  these  techniques  are  extended 
and  applied  to  the  problem  of  forced  vibration  of  a  system  which  has  a 
small  stochastic  parameter.  As  an  example,  the  techniques  are  applied 
to  a  string  of  randcai  density. , 
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1.  Introduction 


The  class  of  problems  involving  stochastic  systems  can  be  classified 
broadly  into  three  categories 

i )  deterministic  System  -  Stochastic  Input 

ii)  Stochastic  System  -  Deterministic  Input 

iii)  Stochastic  System  -  Stochastic  Input. 

Of  these  three  classes,  the  most  highly  developed  analytical  techniques 
are  available  for  i.  A  discussion  of  this  topic  was  recently  given  by 
Ames  [1].  Most  of  the  techniques  for  ii  and  iii  are  approximate,  and  are 
valid  for  systems  with  small  stochastic  parameters.  The  methods  described 
herein  are  applicable  to  the  eigenvalue  problem,  and  the  inhomogeneous 
boundary  value  problem  for  weakly  stochastic  systems.  Although  this  may  seem 
restrictive,  it  does  nevertheless  encompass  a  large  spectrum  of  meaningful 
problems  such  as  free  vibration  or  buckling  or  a  structure  whose  average 
properties  are  known  (such  as  mass  produced  structures ) ,  vibration  of 
stochastic  structures  forced  by  either  random  or  deterministic  loading, 
static  response  of  stochastic  structures  to  either  random  or  deterministic 
loading,  stochastic  control  systems  and  many  others. 

2.  The  Eigenvalue  Problem 

Instead  of  discussing  the  general  problem  immediately,  let  us  consider 
the  following  problem.  We  have  a  taut  string  of  length  L,  and  tension  T. 

The  density  of  this  string  is 

p(x)  =  pQ  [1  +  c(x)J 

where  e(x)  is  a  random  variable  which  is  bounded  by  a  small  number. 

PL | e(x) | >aj  =  0 
and  a  «  1. 
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The  equation  which  governs  the  small  amplitude  free  motion  of  this  string 


2  2 
T  4  -  P  LI  +  c(x>]  14  =  0 

3x2  0  at2 

where  y(x,t)  is  the  displacement  of  the  string  from  its  neutral  position. 

An  eigenvalue  problem  arises  when  y{x,t)  =  Y(x)e1Wtt  and  y(Q,t)  =  y(L,t)  =  0. 


Namely, 


+  X(.l  +  e(x)  J  Y  =  0 


where 


A  .  £ 


We  restrict  our  attention  to  the  fundamental  mode  (although  the  same 
technique  can  be  applied  to  any  mode),  and  realize  that  for  each  function 
e(x)  there  is  a  corresponding  X.  That  is,  X  is  a  functional  of  e(x)  on 
the  interval  [o,Lj.  It  can  be  shown  that  X  is  a  continuous  functional, 


and  therefore  can  be  represented  as  follows 


,  n)e(£)e(n )d(dn 


X=V  \  ^(OetOdC  ♦  \  1  K  (C1n)c(C)e(n)cUdr 

fL  fL 
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Cz  .  f .  Volterra  [2j). 

We  recall  that  |e(x)j  is  uniformly  bounded  by  a  small  number  and,  as 
an  approximation,  truncate  the  above  expression  after  the  first  integral. 


X  =  X  ♦ 
o 
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(£)  is  evaluated  by 


A  is  the  eigenvalue  corresponding  to  e(x)  =  0,  K. 
o  i 

taking  the  functional  derivative  (c.  f.  £2]).  For  this  problem 


v  (r\  -  2lr  qi  2 

K1(U  -  -  X  1“ 

if 


Hence , 


.  n2  2  it2  T  .2 

-7  ~  1  sin  ~ 


e(4)dC  • 


Since  e(f, )  is  a  random  variable,  so  is  A.  The  problem  is  completed  when 
the  statistical  properties  of  A  are  related  to  those  of  e(t, ).  For  the 
remainder  of  this  report  the  expected  value  of  a  random  variable  u  is 
denoted  by<u>  . 


It  follows  immediately  from  (3)  that 
<A7  =  n2/u2  -  ~~  (  S. 


f  Sin2  J&-  <e(0>  d«  . 


In  the  particular  case  tiiat 


^c(C)>  *  o,  we  have 

<A>  =n2/.2 


(A-  <A>  )2  *  f  C  Sin2  j^Sin2  J2.  {c(C)e(n>^  d£dn  .  (5) 

L  ^0 

The  quantity  ^e(C)t(n)^  is  the  autocorrelation  of  e(x)  tind  so  we  related 
the  mean  value  and  standard  deviation  of  A  to  the  mean  value  and  autocorre¬ 
lation  of  r.(x). 

This  illustrative  problem  suggests  the  following  general  approach. 


Suppose 


Mlc(x) J  y ( a ;  +  A  y(x)  =  0 
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where  MLc(x)J  is  a  linear  differential  operator  which  involves  the  small 
random  variable  e(x).  The  boundary  conditions  are  homogeneous  at  x  =  0 
and  x  ■  L.  We  further  suppose  the  eigenvalue  problem  has  solutions  in 
the  neighborhood  of  e(x)  =  0.  Denote  the  smallest  eigenvalue  A,  and 
the  corresponding  value  for  e(x)  *  0.A Q.  Then  by  the  same  reasoning 
as  before  we  can  approximate  A  as  follows 


A  -  A_  + 


KU)  $(OdC. 


from  this  it  follows  that 


< ^  =  Ao  + 


K<0  e(£)  d£ 


and  if  *  0 


(A-  <A>  Y 


rr 
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K(£)  K(n)  ^e(4)£(n)>  d£dn. 


The  successful  application  of  this  method  requires  calculation  of  K(£)» 
but  this  is  not  too  difficult  [2]. 

3.  Inhomogeneous  Boundary  Value  Problem 

Consider  the  following  class  of  problems 

Mle(x)J  y(x)  =  f(x)  (8) 

where  M[e(x)j  is  a  linear  differential  operator  depending  on  the  uniformly 
small  function  c(x).  Furthermore  we  assume  homogeneous  boundary  conditions 
at  x  =  0,  x  =  L.  Such  problems  are  easily  solved  if  one  knows  the  Green's 
function  for  the  problem.  The  Green's  function  is,  however,  dependent  on 
the  function  e(x).  Let  Ge(x,£)  be  the  Green's  function  corresponding  to 
a  particular  c(x).  Then  the  solution  to  (1)  for  that  c(x)  is 


■  r 


G£(x,C)  f(S)  dC. 


out  G^(XfO  is  a  functional  of  e(n)»  and  in  the  case  of  uniformly  small 
c(n)  can  be  approximated 

rL 

Ge<x»£)  i  Gq(x,0  +  \  K(x,£,n)  e(n)  dn  . 

Hence  the  approximate  solution  to  (1)  is 


■j: 


Gq(x,0  f(Od£ 


rr 

^  0  o 


K<x,C,n)  e(n)  f <£)  dndf.  • 


Where  Gq(x,S)  is  the  Green's  function  corresponding  to  e(n)  =  0, 
and  K(x,£,n)  is  the  functional  derivative  of  the  Green’s  function  with 
respect  to  c(n)  evaluated  at  e(n)  =  0.  The  above  formulation  is  applicable 
for  random  f(fc)  as  well,  and  taking  averages  we  obtain 

<y(x)>  =1  Gq(x,0  ^  f(£)^  dc 

J  0 

♦  \  I  K(x,C,n)  ^e(n) f(U>  dndc  (11) 

»  o  •  0 


and,  in  case  ^f(C)>  =  f(0  and  ^e(n)^  *  0 


we  obtain 


(y(x)  =  ^  y(x)^  )4 


I  I  1  I  ^U,s,t)  f(0  t(s)  <c(n)«(t)>  d4dsdnd^  (1/ 
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As  an  example,  consider  the  following  problem 


2 

♦  [l+e(x)]y  =  6(x-L/2) 

dx  (13) 

y(0)  =  y(L)  =  0. 

Physically,  this  corresponds  to  a  harmonic  point  load  applied  at  the 
center  of  a  taut  string  whose  density  is  constant  plus  a  small  random 
variation*  We  further  assume  that 


^c(x)>  =  0 

(14) 

and 

^e(x)  e(y)^  =  Ae  b^x  yL 

(15) 

For  equation  (13)  and  the  boundary  conditions,  one  obtains 


r 


K(x,£,n) 


K 


SinX2|^a~°  LCot(L-0-Cot  Li 

«“«• 

Si~SiTnr<-L~~  CCot(L-x)-CotH 


x<t<n 


x<n<t 

n<x<£ 

(lb) 


t<x<n 


Sint  S in(L-x) 
2STn£ 


Cot  L 


t<n<x 


From  eqs.  (11),  (14),  (12),  (15)  we  obtain 
^  y(x)>  =  0, 


b 


tr  | 
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